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ABSTRACT 



We report the detection of a massive (M^ M >5xlO Mq) molecular/dusty disk of 1.4kpc-size fueling the central engine of the Compact 
Symmetric Object (CSO) 4C 31.04 based on high-resolution (0.5"-1.2") observations done with the IRAM Plateau de Bure interfer- 
ometer (PdBI). These observations allow for the first time to detect and map the continuum emission from dust at 218 GHz in the disk 
of a CSO. The case for a massive disk is confirmed by the detection of strong HCO + (1-0) line emission and absorption. The molecular 
gas mass of 4C31.04 is in the range 0.5xl0 10 -5xl0 l0 M o . While the distribution and kinematics of the gas correspond roughly to 
those of a rotating disk, we find evidence of distortions and non-circular motions suggesting that the disk is not in a dynamically 
relaxed state. We discuss the implications of these results for the general understanding of the evolution of radio galaxies. 

Key words. Galaxies:Individual: 4C 31.04-Galaxies: ISM-Galaxies: kinematics and dynamics-Galaxies: active-Galaxies:nuclei. 



1. Introduction 

Contrary to the well developed radio jets of classical radio galax- 
ies (on scales of ~a few 10-100 kpc), the distinguishing property 
of Compact Symmetric Objects (CSOs) is the small separation 
between their radio-lobes (~a few 100pc-l kpc). Two scenarios 
were proposed to account for the radio structure of CSOs: the 
youth scenario, which contemplates CSOs as young radio galax- 
ies with jets in the making (e.g. Fanti et al. 119951 Owsianik & 
Conway 1998 ), and the frustration scenario, which interpreted 
the absence of large jets as a consequence of the confinement by 
a highly dense massive cocoon surrounding the central engine 
(e.g. O' Dea et al. [19911 Bicknell et al. [19971 Of particular note, 
there is mounting evidence that interaction between the radio 
plasma and the ISM present around the central engines of radio- 
loud galaxies is at work also in CSOs (Morganti et al. 2004 ) . 

A study of the content, distribution and kinematics of cold 
gas via molecular and atomic lines is key to elucidating the 
nature of CSOs and their evolutionary link to other radio- 
loud galaxies. 4C 31.04 is a nearby CSO at z~0.06 at the nu- 
cleus of the giant elliptical MCG 5-8-18, located at the cen- 
ter of a cluster. HI absorption was detected by Van Gorkom et 
al. (119891 ) and Mirabel (1990), and mapped at high-resolution 
by Conway ([1996). The HST color maps obtained by Perlman 
et al. (T2001l )(Pe01) unveil that the optical image of 4C 31.04 is 
permeated with obscuration features. These consist of an elon- 
gated ~l"-central disk oriented along PA~5°-10°, surrounded 
by a spiral-like structure which stretches from the disk over 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



the central ~4" of the image. The central disk is roughly per- 
pendicular to the axis of the double-lobe radio source shown 
in the 5 GHz VLBA image of Giovannini et al. (l200TT) (Gi01). 
This radio source consists of two lobes that are separated by 
~0.07"(80pc) and extend close to the plane of the sky. 

We have probed the gas and dust content of 4C 3 1 .04, by ob- 
serving the HCO + (1-0) and I2 CO(2-l) molecular lines and their 
underlying continuum emission with unprecedented resolution 
(<1") using the IRAM Plateau de Bure Interferometer (PdBI). 
In this Letter we discuss the results obtained from the analysis 
of the PdBI maps of continuum emission detected at 1mm and 
3mm. We also analyze the HCO + (1-0) map of the galaxy show- 
ing emission and absorption of this line. The CO(2-l) molecu- 
lar line has only been detected towards the central offset of the 
galaxy due to insufficient sensitivity. The results obtained show 
the existence of a highly disturbed massive (M ?aj >5xl0 9 MQ) 
molecular disk fueling the central engine of 4C 3 1 .04. 

2. Observations 

We observed 4C 31.04 simultaneously in the HCO + (1-0) and 
I2 CO(2-l) transitions with the PdBI in the A+ configuration 
on February 2006. We used one receiver setting centered on the 
HCO + transition (redshifted to 84.204 GHz) and two overlap- 
ping settings (centered each at 217.433 GHz and 217.873 GHz) 
to sample the 12 CO(2-l) line (redshifted to 217.653 GHz). 
These receiver settings translate into a velocity coverage of 
^OOOkms- 1 at 84GHz and -HOOkms- 1 at 218 GHz. The 
phase tracking center was set at (a^ooo, £2()oo)=(01 h 19 m 35.001 s , 
32°10'50.06") coincident with the position of the radiocon- 
tinuum source core (GiOl). The line visibilities at 84 GHz 
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Fig.l. a) The lmm-continuum PdBI map (at 218 GHz) of 
4C 31.04 (contours: 3<x, 4cr, 7<x to 55<x in steps of 6cr; 
lcr=0.8 mJy beam 1 ) is overlaid on the R-H color image 
of PeOl from HST (color scale). (Aa, A£)-offsets in arc- 
sec are relative to the location of the AGN. See discussion 
in Sect. O for details, b) The HCO + (1-0) line maps of 
4C 3 1 .04 obtained after subtraction of the continuum source 
for the three velocity channels defined in Fig. [3] I(emission)=[- 
300,25] km s _1 (black thin contours; with maximum at B), 
II(absorption)=[25,275]kms _1 (grey contours; with maximum 
at the AGN), III(emission)=[275,900]kms" 1 (black thick con- 
tours; with maximum at A). Levels are 3.5<x, 4cr, and 4.5cr for 
channel I (lcr(I)=0.06Jykms _1 ), -40<x, to -lOOcr in steps of - 
10cr for channel II (lcKII^O.OSJykms- 1 ) and 5cr to 13cr in 
steps of 2cr for channel III ( 1 cr(III)=0.08 Jy km s" 1 ). HCO + con- 
tours are overlaid on the F702W HST image of PeOl (color scale 
in mag.pixeL 1 ). Ellipses represent the beams at 218 GHz(a)) and 
84GHz(b)). See discussion in Sect. l3.2l for details. 



were self-calibrated on the point-like radio-continuum source 
in the 84.402-84.475 GHz band. The l<x-noise for an on- 
source integration time of 9hrs is 0.6mJybeam _1 in 7.5 MHz- 
channels. At 218 GHz the visibilities were calibrated using the 
antenna based scheme. The point-source sensitivity at 218 GHz 



is 1.9 mJy beam 1 in 15 MHz-channels. Maps were obtained 
with uniform weighting and yielded synthesized beams of 
0.50"x0.24" @PA=12° at 218 GHz and 1.20"x0.64" @PA=26° 
at 84 GHz. We have obtained a map of the continuum emis- 
sion at 218 GHz in 4C 31.04 using the 1 GHz- wide image 
side band of the 1mm receivers. Inside this window, there 
is no contribution from line emission. Similarly, a 3mm 
continuum map (at 84 GHz) has been built using velocity 
channels free of line emission in the signal side band of 
the PdBI receivers. The corresponding 1-cr sensitivities are 
-0.8 mJy beam^ 1 and -0.5 mJy beairT 1 at 218 GHz and 84 GHz, 
respectively. Although re-determined in this work, the veloc- 
ity scale is initially referred to the redshift z=0.0592 de- 
rived by Sargent (1973). Luminosity and angular distances are 
D L =260Mpc and D A =230Mpc; the latter gives 1"=1.1 kpc. 

3. Results 

3.1. Continuum maps: detection of a dusty disk 

The 1mm map shown in Fig.QJ reveals a strong emission source 
in the nucleus of 4C 31.04. The source can be decomposed into 
an unresolved central component, and a disk. The disk extends 
over ~1.6"(~1.8kpc) along PA~5°-10°. The measured lmm- 
continuum visibilities can be fitted by a ~40mJy point source 
centered at the AGN, and an elliptical Gaussian source which has 
a deconvolved FWHM-size of l"x0.1" (1.1 kpc x 100 pc) and an 
integrated flux of ~20mJy. The continuum emission at 3mm 
consists of a point source of ~160mJy with no extended disk 
detected to the limits of our spatial resolution and sensitivity. 
The point source detected at 1mm and 3mm is the higher fre- 
quency counterpart of the non-thermal double-lobe radio source 
mapped at cm-wavelengths in this CSO (e.g., GiOl). The spec- 
tral index (a) of the total continuum emission between 84 GHz 
and 218 GHz (with S v ~ v~ a ), confirms the dominance of the 
non-thermal central component: a> 1 . 

The detected 1mm disk is tightly linked to the dusty disk 
identified in the R-H color image of PeOl (see Fig.QJ). The 
dusty disk seems to be mostly edge-on and oriented perpendicu- 
lar to the axis of the double-lobe radio source. However, thermal 
free-free emission from an embedded star forming episode could 
also partly contribute to the disk emission at 1mm. While there 
is no evidence of such a disk at cm wavelengths, it should be 
noted that the published VLBA/VLBI maps of 4C 31.04 would 
not be able to recover any emission on the typical scales of the 
disk. However the non-detection of an extended disk at 3mm ar- 
gues against a significant free-free emission component at 1mm. 
Furthermore, there is no detection of a thermal outflow-like com- 
ponent either at 1mm or at 3mm. 

If the 1mm flux of the disk can be attributed to dust emis- 
sion, a mass budget (Mj„„) can be derived provided that val- 
ues for the temperature (Tdust) an d f° r the emissivity of dust 
(k) are assumed. IRAS observations reveal significant amounts 
of warm dust in 4C 31.04. From the flux ratio measured by 
IRAS S v (100//)/S v (60//)~3.5, we derive a range for T dust =25- 
30 K for a corresponding range of /3-2-1 in the wavelength de- 
pendency of k ~ A~~P . Taking /<-(850yu)=0.9cm 2 gr _1 from Klaas 
et al. ( 120011 ), and T rfHJ ,=30K (i.e., the solution for j3-l) we de- 
rive Md ust ~6xlO 7 M0 for the warm dust. To better estimate the 
total dust mass we have fitted the fluxes measured by IRAS (at 
100// and 60//) and the PdBI (at 1.3mm) using a two compo- 
nent model for the dust temperatures. Any plausible fit to the ob- 
servations gives a value for Mj„„ of (4-7)xlO 8 M0. The higher 
(lower) value of Ma ust corresponds to the fit with y3= 1 .5 (J3-1) 
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Fig. 2. The kinematics of molecular gas along the strip connect- 
ing points A and B (emission peaks for HCO + (1-0)) are illus- 
trated by a position-velocity diagram. Contour levels represent 
emission of HCO + (1-0): 2cr, 2.5<x, 3cr to 7cr in steps of lcr; 
lcr=0.6mJybeam _1 . Color scale (mJybeairT 1 ) represents ab- 
sorption in HCO + (1-0) near the AGN. X-offsets are in arcsec 
relative to the position closest to the AGN. Velocity scale in Y 
axis referred to v o (z=0.0592). 



and temperatures of ~16K(21 K) and ~32K(39K) for the cold 
and warm dust components, respectively. The cold dust compo- 
nent seems to contain the bulk of the mass: M™'^ >0.95xM^'). 

Taking a molecular gas-to-dust ratio ~100 (Seaquist et 
al . 120041 1, we foresee that the total molecular gas mass of the disk 
might be M^ aj ~(4-7)xl0 10 Mo. If we allow for a minor contri- 
bution of free-free emission to the 1mm flux, this can be taken 
as an upper limit to M gas . As argued in Sect. 13.21 the case for a 
massive gas reservoir in 4C 31.04 is nevertheless supported by 
the detection of HCO + (1-0) emission from a rotating molecular 
disk. 



3.2. Line emission: detection of a molecular disk 
3.2.1. Distribution and kinematics 

The channel maps obtained after subtraction of the 3mm contin- 
uum source show that HCO + (1-0) is detected both in emission 
and absorption in 4C 31.04. As shown in Fig.QJ), HCO + (1-0) 
is detected in emission in two velocity ranges: from -300 to 
25kms _1 , and from 275 to 900 kms 1 (channels I and III, re- 
spectively; see Fig. [3). The strongest emission (III) peaks NE 
with respect to the AGN (offset A=[0.1", 0.4"] in Fig.QJi). 
Compared to III, emission in I is >3 times weaker, and peaks 
SW with respect to the AGN (offset B=[-0.6", -0.6"] in 
Fig-lb)- In contrast, the line is detected in absorption from 25 
to 275kms~ 1 (channel II; see Fig. [3}. The HCO + absorption is 
detected as a point source centered at the AGN. These results 
indicate that the distribution of molecular gas traced by HCO + 
has been spatially resolved by the PdBI. 

At first order, the emission and absorption of HCO + seem 
to arise from a rotating disk of ~1.4kpc-size. This is illus- 
trated in Fig. [2] which shows the position-velocity (p-v) plot 
along the strip connecting A and B, i.e., close to the apparent 
kinematic major axis of the disk. By assuming that the total 
width (FWZP) of the HCO + emission profile is twice the pro- 
jected rotational velocity (v rol ) we derive v rot ~500kms _1 for 
an edge-on disk geometry and a value for the systemic veloc- 
ity (v svs ) which is +300kms~ 1 -redshifted with respect to the 




Fig. 3. a) Spectrum of the HCO + (1-0) line emission and absorp- 
tion at the central offset (AGN) of 4C 31.04. A strong absorp- 
tion line is detected against the 3mm continuum source. Velocity 
scale is referred to v , which corresponds to z=0.0592. b) Same 
as a) but for the CO(2-l) line, c) A zoomed view on the spec- 
trum shown in a) helps to delimit intervals I-to-III and derive a 
new value for v 5US . 



redshift initially assumed in this work (z=0.0592). The HCO + - 
based redshift (z(HCO + )=0.0602+0.0002) agrees with the val- 
ues recently derived from optical spectroscopy observations of 
4C 31.04 (z=0.0600+0.0001; Marcha et al . 1X9951 . 

If circular motions completely ruled the gas kinematics, the 
absorption profile against the AGN (this being the dynamical 
center) should be centered at v jyj . This is however at odds with 
observations: as shown in Fig. [5] the gas seen in absorption in 
channel II is ~150kms _1 blue-shifted on average with respect 
to z=0.0602 (similarly to the HI absorption of Conway 1996). 
This result can only be reconciled if the gas causing the ab- 
sorption is outflowing from the nucleus or subject to strong non- 
circular motions, thus partly invalidating the initial premises. An 
alternative solution adopting the center of the absorption pro- 
file as v svs would only enhance the asymmetry in the HCO + 
emission profile. Furthermore, there is evidence of distortions 
in the molecular disk from the observed distribution. The SW 
side of the molecular disk is offset from the dust disk: the HCO + 
disk is tilted towards a larger PA (Fig.QJ)). A conspiracy of ab- 
sorption and emission may explain the observed anticorrelation. 
However, independent evidence of distortions is also found on 
the opposite side of the disk: whereas the bulk of the HCO + 
emission NE follows the dust disk along PA=5°-10°, the molec- 
ular disk is seen to thicken along its minor axis in this region 
(Fig.[D>). 

While the reasons explaining the observed distortions in the 
distribution and kinematics of the gas cannot be elucidated with 
these observations, these results further illustrate that the rotating 
molecular disk of 4C 3 1 .04 is not dynamically relaxed. 
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3.2.2. Masses 

The HCO + line in emission probes the dense phase 
(n(H 2 )>10 4 cnr 3 ) of the molecular disk of 4C 31.04. Without 
further constraints at hand, an estimate of the mass of dense 
molecular gas (Md mse ) can be derived from the luminosity of the 
line (L[HCO + ]=4.2xl0 8 Kkms~ 1 pc 2 ) if we assume low opaci- 
ties, LTE conditions, an abundance for HCO + (X[HCO + ]), and 
an excitation temperature (T fA ) for the line. For T„=10K and 
X[HCO + ] = 10~ 8 , we derive M dense ~5xlO 8 M . This is probably 
a lower limit, because the HCO + line is likely optically thick (see 
Fig. 0). Alternatively, we can derive Mdense via a L(HCO + )-to- 
Mdense conversion factor if we assume that the emission comes 
from a clumpy medium of optically thick virialized clouds. The 
value of the conversion factor would depend on the physical pa- 
rameters of the clouds through the ratio ~n^ 2 (H2)/Ti,(HCO + ), 
where n(H2) is the molecular gas density and T/,(HCO + ) is the 
emergent line brightness temperature of the clouds. Through 
LVG simulations, assuming a range for T^=20-to-60K and 
n (H2)=10 4 -to-10 5 crrr 3 , we derive an average conversion fac- 
tor M d( ,„ je /L(HCO + )~10M o (Kkms- 1 pc 2 )- 1 . This would im- 
ply that M tto „ ( ,~4xl O 9 M . Taking HCO + (1-0)/CO(1-0)=0.05 
(Gracia-Carpio et al. 120061 1 we derive Mg flJ /M^„ se ~10, and thus 
estimate M gas to range between ~5xl0 9 M G and ~4xlO lo M . 

The sensitivity reached in CO(2-l) was insufficient to map 
the emission of this line in 4C 31.04, therefore M gas cannot 
be derived from these data. The line is detected in absorption 
and emission only towards the AGN (Fig. [3]). However, Ocana- 
Flaquer et al.(in prep) recently detected in 4C 31.04 the emis- 
sion and absorption of CO(1-0) and CO(2-l) using the IRAM 
30m telescope, improving previous searches for CO emission 
(O'Dea et al. 120051 Evans et al. 120051 Using a CO-to-H 2 con- 
version factor N(H 2 )/Ico=2.2xl0 20 cm _2 K _1 km -1 s (Solomon 
& Barrett[199T), Ocana-Flaquer et al. estimate M^ fls ~5xlO 9 M , 
a value within the range derived from HCO + . Note that due to the 
superposition of absorption and emission in the profiles, the val- 
ues derived from molecular line emission underestimate M gas . 

Taken the estimates from dust continuum and line emis- 
sion together, we conclude that M gas lies in the range (0.5- 
5)xl0 10 M Q . This is much larger than the mass of neutral hy- 
drogen derived by Conway ( 1996 ) from HI line data (~1O 8 M ) 

3.3. Line absorption 

Assuming LTE conditions we can derive the column densities 
of HCO + (Nhco+) and CO(N C o) from the the J->J+1 absorption 
lines detected towards the center of 4C 31.04 according to: 

8ttv 3 Q(T ex )exp(Ej/kT ex ) f 

N mo l = -= I TdV 

c s A u ig u (1 -exp(-hv/kT ex )) J 

where Ey is the energy of the lower level of the transition, T ex is 
the excitation temperature, Q(T ex ) = 2y/(27/+l) exp(-Ejj/kT ex ) 
the partition function and A,,/ the Einstein coefficient of the line. 
A lower limit to the velocity integrated optical depth of the line 
J TdV, and thus to N mo i, can be derived assuming that the cov- 
ering factor / of the absorbing molecular gas with respect to 
the extent of the radio source in 4C 3 1 .04 is unity. In this case 
we replace J TdV by T obs Av, where r oij =ln(T^/T^) with T c mb 
and Tfj. being the main brightness temperature of the contin- 
uum and the line, respectively, and Av the Gaussian width of 
the line. From r & s we estimate a lower limit to / which is 
similar for the two lines (/ >l-e~ T "'" ~0.27), suggestive of a 
common value for /. Assuming f-1 and 7" ex =10K for the two 



lines, we derive N(HCO + )/N(CO)=1.6xlO- 3 . While this value 
is surprisingly high relative to galactic standards (see Table 3 of 
Wiklind & Combes 119951 [Wi95]), it likely overestimates the 
true HCO + /CO abundance ratio for two reasons. First, T ex is 
probably smaller for the HCO + line than for the CO line. Second, 
the chances for the CO absorption to be underestimated due to 
corruption by emission within the beam are much higher than 
for the HCO + absorption (see discussion in Wi95). 

4. Conclusions 

The high-resolution and high-sensitivity capabilities of the PdBI 
have allowed for the first time to detect and map a massive 
molecular/dusty disk fueling the central engine of a CSO like 
4C 31.04. Based on the estimates derived from the contin- 
uum and line emissions, the molecular disk mass M gas lies in 
the range (O.5-5)xlO K) M . Even if we adopt the lower bound 
to M gas , these observations reveal the existence of significant 
amounts of molecular gas in the central 1 .4 kpc disk of 4C 3 1 .04, 
comparable to those identified in ULIRGs (~10 9 -5xl0 10 M o ; 
Sanders et al. 1991), and much larger than the masses of cold 
gas detected in normal elliptical galaxies (~10 7 -10 8 M o ; Lees et 
al. 1 1 99 11 >. Of particular note, the large value of M gas derived for 
4C 31.04 would be roughly consistent with the scenario of con- 
finement in this radio source. Using simple analytical models, 
Carvalho ( 1998 ) predicts that a clumpy medium with M gflJ ~ 10 9 - 
1O 1O M can confine a radio source on the scale of 0.5-1 kpc. In 
the case of 4C 3 1 .04, the bulk of M gas is seen to be in a disk 
roughly perpendicular to the radio source axis. This particular 
geometry would decrease the efficiency of frustration, however. 

Within the frustration scheme, it is foreseen that a disturbed 
ISM around the radio AGN would reflect the interaction between 
the radio plasma and the dense confining cocoon. In the case 
of 4C 3 1 .04, there is evidence that a radio lobe-ISM interaction 
may be at work. The HST image of PeOl shows cone-like fea- 
tures aligned with the radio axis of the source (see Fig. QJ)), 
the likely signature of gas shocked by the jet (e.g., Labiano et 
al. 120031) . Moreover, the distribution and kinematics of molecu- 
lar gas probed by the HCO + observations presented in this work 
illustrate that the detected rotating disk is not in a fully relaxed 
state. The revealed distortions may reflect that the disk is still set- 
tling after a merger or an event of gas accretion. Alternatively, 
the jet and the cone-like features may be interacting with the disk 
and thus produce the reported distortions. 

Obtaining a better constrain to M gas in 4C 31.04, and in 
CSOs in general, is key to exploring the nature of these sources. 
As illustrated in this work, the central gas disks of radio-galaxies 
are likely dominated by the molecular phase. In order to take 
full advantage of molecular line data, high spatial resolution and 
high sensitivity observations are needed to reduce at most the 
bias due to the confusion of absorption and emission in the beam. 
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